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SUMMARY 

This  investigation  was  performed  to  determine  the  feasibility  of 
obtaining  lateral  frequency  response  and  stability  derivatives  utilizing 
the   forced  oscillation  technique.     Tests  were   conducted  in  a  North 
American  Navion  at  one  altitude,   airspeed  and  gross  weight  condition,, 

Both  rudder   and  ailerons   were   oscillated  sinusoidally  with  a 
simple  mechanical  oscillator  and  aircraft  responses  were  telemetered 
and  recorded  on  magnetic  tape     A  technique  was  developed  for  elimin- 
ating distortion  of  the  input  sine  wave  by  utilizing  a  low  frequency  band 
pass  filter  in  the  reduction  of  flight  data.     Frequency  responses  were 
determined  for       Cf)     ,       ~\jf    ,    and      /?     for  both  the  rudder   and  the 
aileron  oscillations;   and  these  were  used  to  determine  the  aircraft's 
stability  derivatives.     Frequency  responses  were  calculated  from  the 
experimentally  determined  derivatives  and  were  compared  with  the 
actual  frequency  responses  to  illustrate  the  accuracy  of  the  method 
employed  to  extract  the  stability  derivatives. 

The  methods  employed  proved  entirely  feasible  for  this  aircraft 
at  the  single  flight   condition  and  apparently  were  applicable  to  many 
current  production  aircraft.     The   frequency  responses  showed  little 
random  error,   and  the  responses  calculated  from  the  experimentally 
determined  derivatives  closely  matched  the  actual  responses. 
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INTRODUCTION 

The   design  of  modern  high   speed  air cr ait  and  the   automatic 
control  systems  with  which  nearly  all  are  equipped  requires   an 
accurate  evaluation  of  their  response  characteristics.     In  the   early 
design  stages  the  only  method  available  to  the  designer  involves   the 
estimation  of  the  aircraft's  stability  derivatives  and  the  solution  of 
the  assumed  equations  of  motion.     Because  of  the  inherent  difficulties 
involved  in  the  analytic  determination  of  stability  derivatives,   large 
errors  may  be  present  in  the  derived  response  characteristics. 
Dynamic  flight  testing  of  the  prototype  aircraft  to  accurately  determine 
its  stability  derivatives  and  the  coefficients  of  its  transfer  function  is 
becoming  increasingly  important. 

Dynamic  flight  test  data  is  usually  obtained  by  one  of  two 
methods.     The  most  popular  of  these  involves  the  determination  of  the 
transient  response  of  the  aircraft  to  known  control  inputs.     While  the 
data  thus  obtained  is,   in  itself,  not  in  the  best  form  for  analysis,    the 
recent  development  of  efficient  computing  machines  has  made  it  possi- 
ble to  convert  such  a  transient  response  to  frequency  response   form 
where  it  may  be  more  readily  analyzed.     The  second  method  involves 
the  direct  determination  of  the  aircraft's  frequency  response   by  the 
steady  state  forced  oscillation  technique.     While  considerably  more 
flight  time  is  required  in  this  method,   data  reduction  is  greatly  simpli- 


fied ;    and  good  results  may  be  obtained  without  the  use  of  automatic 
computers. 

The  determination  of  an  aircraft's  lateral  frequency  response  by 
the  forced  oscillation  technique  may  be  accomplished  by  oscillating 
either  of  the  lateral  controls  with  a  sinusoidal  forcing  function  while 
keeping  the  other  locked.     From  simultaneous  recordings  of  the   air- 
craft  responses  in      0     ,      Iff     and     /3      and  the  control  inputs   on 
separate  time  traces,   it  is  possible  to  determine  amplitude  ratios  and 
phase  angles  for  each  response  over  a  range  of  oscillation  frequencies 
for  the  steady  state  condition.     From  these  amplitude  and  phase  spectra 
it  is,   in  theory,   possible  to  determine  from  either  the  aileron  or  rudder 
data  all  the  aircraft's  aerodynamic  derivatives,   except  of  course   the 
control  derivatives  for  the  locked  surface,   providing  the  weight   and 
moments  of  inertia  of  the  aircraft  are  known.     In  practice,  however,   it 
is  frequently  difficult  to  determine  all  the  derivatives  accurately  since 
some  of  them  have  little  effect  on  the  frequency  response   and  thus   are 
difficult  to  determine  from  it. 

While  considerable  effort  has  been  expended  on  this  type  of  testing, 
especially  at  the  Cornell  Aeronautical  Laboratory,  the  accurate  deter- 
mination of  derivatives  without  harmonic  analysis  of  the  time  traces  has 
required  very  pure  sinusoidal  control  inputs,,     This  requirement  has  led 
to  complicated  instrumentation  which  has  made  the  method  applicable 
only  to  research  type  testing. 


This  investigation  was  conducted  to  determine  the  applicability 
of  the  method  to  more  simplified  testing  techniques  utilizing   a  simple 
mechanical  oscillator,  telemetering  equipment  for  data  recording,   and 
an  ultra -low  frequency  band-pass  filter  to  eliminate  sine  wave  distor- 
tion.    The  flight  tests  were  conducted  in  a  North  American  Navion  at 
the  Forrestal  Research  Center,   Princeton  University,   during  the 
winter  and  spring  of  1957. 


EQUIPMENT 

Test  Airplane 

The  test  airplane  was  a  North  American  Navion,   which  is  an  all 
metal,  four  place,   low  wing,   single  engine  airplane.     The  Navion  is 
equipped  with  a  Continental  E-185  engine,   rated  at  185  HP  at  2300 
RPM  and  29  inches  of  MP,   at  sea  level.     The  wing  is  unswept,   and  the 
tail  configuration  conventional .     The  control  surfaces  are  all  of  beaded 
skin  construction,  and  have  no  trailing  edge  extensions.     Trim  tabs 
are  located  on  the  elevator,   aileron,   and  rudder. 

The  test  airplane  was  modified  by  the  removal  of  the  two  rear 
cockpit  seats  to  provide  space  for  the  test  instrumentation.     The  test 
airplane  had  a  gross  weight  of  2870  pounds  with  full  service,  pilot  and 
co-pilot,   and  all  test  equipment. 

A  three  view  drawing  of  the  Navion  is  shown  in  Fig.  1  and  the 
general  specifications  of  the  test  airplane  are  tabulated  in  Table  I. 
Fig.  2  is  a  photograph  of  the  test  airplane. 

Power   Supply 

Power  available  for  the  instrumentation  in  the  test  airplane  in- 
cluded both  12  and  28  volt  direct  current  supplied  by  a  28  volt  engine 
driven  generator  and  two  12  volt  batteries  connected  in  series.  A  500 
"VA  voltage  and  frequency  controlled  inverter  was  installed  behind  the 
co-pilot's  seat  to  provide  115  volt,  400  cycle,  three  phase  alternating 
current  power  to  the  gyros  and  telemeter  transmitter. 


Considerable  effort  was  expended  in  obtaining  an  inverter  which 
could  be  frequency  controlled,    since  it  was  very  important  that  the 
three  phase  power  to  the  rate  gyros  be  constant  at  the  calibration 
frequency  of  400  cycles. 

A  circuit  diagram  of  the  aircraft   28  volt  system  and  the  inverter 
system  is  shown  in  Fig.   3. 

Aircraft   Instrumentation 


Standard  aircraft  instruments  were  used  to  determine  indicated 
airspeed  and  altitude.     Measurement  of  the  following  quantities   was   re- 
quired:    sideslip  angle,   yaw  rate,  and  roll  rate  with  respect  to  wind  or 
stability  axes,   and  aileron  and  rudder  position  angles.     In  choosing  the 
type  of  components  needed  to  measure  the  above  quantities,  an  attempt 
was  made  to  utilize  instrumentation  compatible  with  the  telemetering 
equipment.     For  this  reason,  the  transducers  selected  were  of  the 
potentiometer  type,  not  only  because  they  were  capable  of  providing 
signals  of  sufficient  magnitude,  but  also  because  their  excitation  source 
readily  provided  zero  and  full  scale  reference  voltages. 

In  view  of  some  past  history  of  various  amounts  of  coupling  be- 
tween the  response  of  the  aircraft  and  certain  of  the  instrumentation 
components  in  this  type  of  flight  testing,   considerable  attention  was 
given  to  choosing  transducers  with  satisfactory  frequency  response 
characteristics  in  order  to  avoid  correction  of  response  data  for  the 
dynamics  of  the  recording  instruments.     In  general,  it  was  desirable 


that  the  natural  frequencies  of  the  rate  gyros  and  the  sideslip  vane  be 
large  as  compared  to  the  expected  forced  oscillation  frequencies  of  the 
aircraft. 

Sideslip  angle  was  sensed  by  a  light  vane  mounted  on  a  boom 
approximately  one  maximum  fuselage  diameter  ahead  of  the  left  wing 
tip.     The  moment  of  inertia  of  the  vane  was  very  low,   and  no  gearing 
was  introduced  in  the  linkage  in  order  to  keep  the  natural  frequency  and 
damping  within  acceptable  limits.     An  Electro -Mechanical  Laboratory 
ultra-low  torque  potentiometer  with  .008  inch-ounces  of  friction  torque 
was  obtained  for  use  with  the  sideslip  vane. 

In  order  to  reduce    boom  resonance,   which  had  been  experienced 
in  previous  tests  with  the  Navion,   the  boom  was  stiffened  by  three  angle 
strips  equally  spaced  around  the  periphery,   and  running  the  full  length 
of  the  boom.     The  sideslip  vane  and  the  boom  are  shown  in  the  assem- 
bled condition  in  Fig.  4. 

Aileron  and  rudder  angle  were  measured  by  10,  000  ohm  Model 
5301   single  turn  Helipot  precision  potentiometers.,     Due  to  the  nature  of 
the  differential  aileron  system  of  the  Navion,   it  was  not  considered  de- 
sireable  to  pick  off  aileron  angle  at  the  aileron  itself.     Accordingly,   the 
pickoffs  were  made  from  the  control  cable  movement  just  inboard  of 
the  right  aileron.     A  position  close  to  the  aileron  was  chosen  in  order 
to  avoid  the  effects  of  cable  stretchy     The  rudder  deflection  potentiome- 
ter was  mounted  on  the  bottom  of  the  rudder  and  was  actuated  by  a 
cable  under  spring  tension. 


The   measurement  of  yaw  and  roll  rates  was  accomplished  with 
type  3612  Giannini  rate  gyros.     This  type  of  gyro  operates  on  a  power 
supply  of  115  volt,  400  cycle,  three  phase  AC  to  the  motor,   and  115  volt 
single  phase  AC  to  the  heaters.     The  transducer  is  of  the  potentiometer 
type,   with  a  DC  resistance  of  5250  ohms.     These  gyros  were  obtained 
on  loan  from  the  Instrumentation  Laboratory  of  the  U.  S„  Naval  Air 
Test  Center,   NAS  Patuxent  River,   Maryland. 

Oscillating   Mechanism 

Design  of  the  control  driving  mechanism  was  based  on  the  follow- 
ing criteria:    (a)  The  system  must  generate  a  sinusoidal  output  of  con- 
trollable frequency  which  may  be  applied  to  either  the  rudder   or   the 
aileron  control;    (b)  The  driving  link  for  the  driven  control,   and  the 
locking  mechanism  for  the  undriven  control  must  each  include  a  quick 
disconnect,  not  subject  to  binding  when  under  load;    (c)  Driving  and  lock- 
ing links  must  include  a  trim  adjustment;    (d)  Inter -unit  rigidity  and 
close  tolerances  are  required  due  to  the  cyclical  varying  load;    and  (e) 
The  unit  must  be  compatible  with  the  power  and  space  limitations  of  the 
Navion. 

These  conditions  were  met  in  the  following  manner.     The  rotation 
of  a  400  cycle,   1/2  HP  DC  motor  was  converted  to  a  pure  sinusoidal 
translational  motion  by  a  scotch  yoke,   with  suitable  intermediate  gear- 
ing,,    The  driving  rod  to  the  rudder  could  be  disconnected  from  the 
scotch  yoke  plate  and  bolted  to  the  mounting  plate  to  provide  a  positive 
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rudder  lock.     The  rod  from  the  scotch  yoke  plate  was  fitted  inside   a 
tube  from  the  rudder.,     By  use  of  a  large  handled  set  screw  the  rod  and 
tube  could  be  rigidly  locked  at  any  desired  trim  position,   and  quickly 
unlocked  without  possibility  of  binding.     The  aileron  drive  was  obtained 
by  leading  light  aircraft  control  cable  around  pulleys  from  the  horizon- 
tally translating  shaft  of  the  scotch  yoke  to  a  set  of  vertically  tran6lat« 
ing  turnbuckles,   from  which  a  chain  drove  a  sprocket  wheel  mounted 
on  the  control  column.     The  sprocket  was  secured  to  the  control  column 
by  a  large  handled  set  screw  arrangement  which  allowed  aileron  trim 
adjustment  and  a  quick  disconnect.     To  lock  the  ailerons,   the  drive  was 
opened  at  the  turnbuckle  and  a  cable  bolted  to  the  base  plate  substituted 
for  the  driving  cable  . 

The  motor,   gear  box,   and  driving  and  locking  unit  was  built  up  on 
one  large  section  of  1/4"  aluminum  plate,   which  provided  a  very  rigid 
base.     The  complete  mechanism  could  be  installed  or  removed  as  a  unit, 
the  plate  being  bolted  to  extrusions  riveted  to  the  siffeners  under  the 
deck  on  the  starboard  side  of  the  cockpit.     In-flight  shifting  from  the 
aileron  drive  to  the  rudder  drive  involved  about  five  minutes. 

Control  of  the  driving  motor  was  obtained  by  the  use  of  an  ampli- 
dyne  mounted  behind  the  co-pilot.     Through  control  of  the  field  of  the 
amplidyne,   whose  output  provided  the  armature  current  to  the  driving 
motor,  any  speed  from  zero  to  maximum  rated,   including  reversing, 
could  be  obtained.     Indication  of  the  operating  frequency  was  not  directly 


available   in  the  cockpit,  but  could  be  obtained  approximately  from  a 
calibrated  voltmeter  across  the  amplidyne  signal  circuit.     At  the  lower 
frequencies,   the  torque   (proportional  to  the  armature  current)  became 
marginal,   and  was   evidenced  as  a  distortion  to  the   sine   wave  of  the 
driven  control  as  the  aerodynamic  load  varied  with  control  displacement,, 
A  worm  and  gear  was  incorporated  in  the  gear  box  to  obtain  irreversi- 
bility between  the  motor  and  the  control  surfaces,  but  the  variation  of 
the  frictional  load  on  the  worm  was  sufficient  to  cause  this  low  frequency 
distortion.     The  aircraft  instrumentation  is  shown  in  Fig.  6  and  a  close 
up  of  the  oscillating  mechanism  is  shown  in     Fig.  7. 

Signal   Circuit   Design 

The  telemeter   transmitter  was  designed  to  accept  voltages  from 
the  transducers  from  zero  to  five  volts.     It  was  thus  desirable   for  the 
full  scale  voltages  on  the  potentiometers  to  be  maintained  at  five  volts. 
Variations  in  this  full  scale  value  were  permissable,  providing  that  the 
variation  on  all  the  potentiometers  was  the  same,   since   the  telemeter 
corrected  all  voltages  to  a  percentage  of  the  full  scale  reference  voltage. 
The   potentiometer  used  in  the  sideslip  vane  was  a   360   degree    single 
turn  model,  and,  since  only  40  degrees  of  variation  in  sideslip  angle  were 
expected^t  was  necessary  to  apply  45  volts  excitation  to  this  potentiome- 
ter in  order  that  its  output  might  have  a  five  volt  variation.     The  use  of 
a  45  volt  DC  signal  circuit  supply  was  thus  dictated.     It  was  also  desira- 
ble to  maintain  the  signal  voltage  to  the  other  potentiometers  at  five  volts 
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and  introduce  a  biasing  circuit  to  provide  a  variation  of  sideslip  output 
voltage  from  zero  to  five  volts  rather  than  a  five  volt  variation  about 
some  intermediate  voltage. 

The  signal  circuit  design  based  on  the  preceding  requirements  is 
shown  in  Fig.   5.     The  essential  element  of  the  circuit  is  the  bias   sup- 
ply voltage  potentiometer,  pot  "A".     With  the  sideslip  vane  against  the 
stop  in  one  direction,  pot    "A"  was  adjusted  so  that  the  output  voltage 
of  the  sideslip  pot  was  zero.     The  vane  was  then  deflected  to  the  oppo- 
site stop,   and  the  output  voltage  of  the  sideslip  pot  was  recorded. 
This  value  was  slightly  less  than  five  volts.     Pot  "B"  was  then  adjusted 
to  bring  the  full  scale  reference  voltage,   and  thus  also  the  voltage   in 
the  remaining  signal  circuits  in  agreement  with  the  full  scale  sideslip 
voltage.     The  procedure  was  repeated  two  or  three  times  as  necessary, 
in  order  to  nullify  any  loading  on  the  system.     Once  the  final  adjustment 
was  made,  the  zero  to  full  scale  voltage  for  each  quantity  to  be  meas- 
ured was  exactly  the  same,   and  small  variations  in  signal  voltage  pro- 
duced the  same  variation  for  all  circuits.     Pots    "A"    and   "B"   were 
locked  and  required  no  further  adjustment  for  the  remainder   of  the 
project. 

During  the  initial  flight  and  transmission  tests,   considerable  AC 
ripple  was  found  to  enter  the  DC  signal  voltages,   very  likely  from  the 
400  cycle  inverter.     Condensers    "C"  and   "D"  in  Fig.   5  were  installed, 
and  were  found  to  reduce  the  ripple  to  a  negligible  value. 
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Telemeter   System 

The  telemeter  unit  was  of  the  pulse  width  type  and  was  designed 
and  built  by  the  Applied  Science  Corporation  of  Princeton  (ASCOP). 
This  unit  was  the  key  element  in  the  instrumentation  of  this  project. 
The  basic  functional  components  of  this  pulse  width  system  are  trans- 
ducers (all  of  the  potentiometer  type);   a  time  division  multiplexer,   or 
commutator;   an  electronic  coder,   or  keyer  to  put  the  intelligence  into 
pulse  form;   an  RF  transmitter;   and  ground  station  equipment  for   re- 
ceiving,  decoding,   and  recording  the  data  in  usable  form. 

The  airborne  component  of  the  system  is  designated  the  ASCOP 
D  Series  PW  Multicoder,   and  the  ground  station  equipment  the  ASCOP 
M  Series  PW  Ground  Station. 

A  unique  feature  and  very  great  advantage  of  pulse  width  data 
systems  is  that,   due  to  the  multiplexing  or  time  sharing  operation,   all 
data  channels  traverse  the  same  path  during  transmission.     Since  zero 
and  full  scale  references  are  transmitted  with  the  data,   this  results  in 
a  calibration  signal,   effective  for  all  channels,   being  transmitted  20 
times  each  second.     These  calibration  signals  allow  the  removal  of 
almost  all  of  the  system  error,   with  the  exception  of  errors  introduced 
in  the  instrument  pickups  and  the  mechanical  switch,   and  they  make 
possible  the  use  of  this  telemetering  system  in  tests  of  long  duration 
and  varying  environmental  conditions. 
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Magnetic   Tape   Recorder 

An  Ampex  Model  309C  dual  track  three  speed  tape  recorder  was 
an  integral  component  of  the  ground  station  equipment,   and  provided  a 
convenient  means  for  recording,    storage,   and  playback  of  the  standard 
pulse  width  modulated  data.     Provisions  were  also  included  for  placing 
voice  transmissions  directly  on  the  magnetic  tape  from  either  or  both 
the  aircraft  and  the  ground  station. 

Data  Recorder 

Qualitative  analysis  of  the  flight  data  during  the  actual  flights  and 
quantitative  analysis  of  the  recorded  signals  was  made  on  a  Sanborn 
Model  154-100D  four  channel  recorder.     This  instrument  had  a  very 
high  natural  frequency  (42  cps),   and  the  low  frequency  response  was 
flat  to  zero  cps.     Different  paper  speeds  and  variable  gain  were  possible 
and  permitted  expansion  of  the  time  trace  and  amplification  of  the  signal 
to  nearly  any  desired  value. 

Filter 

In  order  to  study  the  response  of  the  aircraft  to  a  true  sinusoidal 
driving  function,   all  forcing  functions  and  the  corresponding  responses 
were  filtered  to  remove  all  but  the  fundamental  frequency.      The   filter 
employed  was  the  model  330-A  KROHN-HITE  Ultra-Low  Frequency 
Band-Pass  Filter.     This  is  an  adjustable  filter  whose  gain  is  unity  in 
the  pass  band  and  drops  outside  the  pass  band  at  the  rate  of  24  db  per 
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octave.     The  filter  combines  a  high  pass  filter  and  a  low  pass  filter, 
each  of  which  may  be  continuously  tuned  from  0.02  to  2,  000   cpe. 
Since  the  filter  will  handle  an  input  signal  of  only  up  to  10  volts   rms 
without  overloading,   a  voltage  divider  was  applied  to  the  telemeter  out- 
put,  which  had  a  full  scale  range  of  up  to  100  volts.     For  this  applica- 
tion the  filter  was  used  as  a  narrow  pass  band  with  the  high  and  low 
dials  set  at  the  same  frequency.     Used  in  this  manner,   the  filter   did 
not  alter  the  relative  amplitude  or  phase  relationships  between  the 
input  and  the  output  at  a  specific  frequency. 

A  photograph  of  the  ground  instrumentation  is  shown  in  Fig.   8. 

Instrument    Calibration 


Six  items  required  calibration  in  the  instrumentation    system. 
These  were  the  two  rate  gyros,  the  two  control  deflection  potentiome- 
ters,  the  sideslip  vane  potentiometer,   and  the  signal  voltage  to  the 
amplidyne  which  was  roughly  a  measure  of  the  frequency  of  control 
oscillation.     The  aircraft  altitude  and  airspeed  indicating  systems  were 
not  calibrated  because  previous  calibrations  indicated  that  the  errors 
here  were  negligible  for  the  conditions  of  this  investigation. 

The  two  rate  gyros  were  calibrated  on  a  pendulum  which  consisted 
of  a   3/4*   diameter  aluminum  shaft  and  a  12"   diameter  instrument  pan. 
The  pendulum,  the  overall  length  of  which  was   121.  75"  was  pinned  to 
the  overhead  truss  structure  in  the  Forrestal  hangar.     A  scale   was 
mounted  with  brackets  on  a  small  table  in  such  a  way  that  displacement 
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of  the  pan  from  its  rest  position  could  be  measured.     The  gyros  were 
mounted  one  at  a  time  on  a  wooden  base  plate  which  was  clamped  to  the 
pendulum  pan,   and  their  potentiometer  outputs  were  connected  to  a 
Brush  Recorder .     The  pendulum  was  displaced  from  its  rest  position 
and  allowed  to  swing  freely.     Since  there  was  an  excessive  amount  of 
friction  in  the  system,   swing  amplitudes  were  measured  on  two  succes- 
sive swings  in  the  same  direction  and  averaged.     The  two  maximum 
voltages  recorded  during  the  swing  were  also  averaged  and  the  maximum 
rate  attained  during  the  swing  was  determined  from  the  following 
relation: 

where   P  was  the  natural  period  of  the  pendulum  whose  measured  value 
from  the  trace  on  the  Brush  recorder  was    3„  33  sec. ,  1  was  the  length 
of  the  pendulum,  and    d    was  the  average  maximum  displacement  of  the 
pan  from  its  rest  position.     The  pendulum  was  allowed  to  swing  until 
once  again  at  rest  with  readings  taken,  as  described  above,   over  the 
entire  range  of  gyro  rates.     The  angular  rate  gradient  of  both  gyros 
was  found  to  be  linear  and  equal  to  .  610  degrees/sec.  per  percent  full 
scale  potentiometer  voltage.     The  gyro  calibration  curves  are  shown  in 
Fig.  9<>     It  was  found  that  the  gyros  were  very  sensitive  to  the  frequency 
of  the  400  cycle  supply  voltage  and  as  a  result  great  care  was  taken 
during  calibration  and  also  during  the  actual  flight  test  program  to 
maintain  the  frequency  of  the  supply  at  exactly  400  cycles  per  second. 


15 


The  control  deflection  potentiometers  were  calibrated  by  measur- 
ing the  resistance  between  the  output  terminals  over  the  full  range  of 
control  deflections.     The  actual  control  deflections  were  measured  with 
a  vernier  protractor  capable  of  accuracy  to  0.1  degree     The  rudder  de- 
flection gradient  was  found  to  be  linear  and  equal  to  .405  degrees  per 
percent  full  scale  potentiometer  voltage.     The   Navion  aircraft  is 
equipped  with  differential  ailerons  and  it  was  expected  that  the  aileron 
calibration  would  therefore  be  non-linear.     However,   using  average 
aileron  angle,  it  was  found  that  over  the  small  deflection  range  of  +10 
to    -10  degrees  desired  during  the  flight  test  program  the  aileron  deflec- 
tion gradient  was  linear  and  equal  to  .1933  degrees  per  percent  full 
scale  potentiometer  voltage. 

The  sideslip  vane  potentiometer  was  calibrated  in  the  same  man- 
ner as  the  control  deflection  potentiometers,   using  a  protractor   and 
measuring  the  resistance  between  the  output  terminals  with  an  ohm- 
meter.     The  calibration  curve  is  shown  in  Fig.  10.     The  sideslip 
gradient  was  found  to  be  linear  and  equal  to  .468  degrees  per  percent 
full  scale  potentiometer  voltage. 

The  telemeter  equipment  corrected  its  output  to  percent  full  scale 

voltage  and  since  full  scale  voltage  on  all  the  potentiometers  varied  by 

the  same  increment  as  the  reference  voltage  for  small  fluctuations  of 

the  latter  value,   it  was  unnecessary  to  calibrate  the  potentiometers 

through  the  telemeter  equipment. 

The  speed  control  on  the  oscillating  device  was  calibrated  roughly 

in  flight  by  timing  the  frequency  of  control  oscillations  and  recording 
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the  voltmeter  reading  across  the  control  field  of  the  amplidyne  genera- 
tor.    This  speed  control  was  used  only  to  obtain  order  of  magnitude 
settings,   the  actual  frequency  of  any  particular  oscillation  being  deter- 
mined from  the  trace   on  the  Sanborn  Recorder. 

Careful  consideration  was  given  to  the  possible  errors  which 
might  be  introduced  in  the  data  by  the  dynamics  of  the  recording  instru- 
ments.    The  natural  frequency  of  the  rate  gyros  was   17  cycles  per 
second  and  their  damping  ratio  was  between  .45  and  .7.     The  natural 
frequency  and  damping  ratio  for  the  sideslip  vane  were  experimentally 
determined  in  flight  to  be  15.44  and  0.  081  respectively.     Since  the 
natural  frequency  of  the  aircraft  was  approximately  0.4  cycles   per 
second  and  the  maximum  forced  oscillation  frequency  was  not  expected 
to  exceed  one  cycle  per  second,  the  attenuation  of  the  signals  due  to  the 
dynamics  of  the  instruments  was  effectively  zero.     While  a  small  error 
in  phase  angle  measurement  was  introduced  at  the  higher  frequencies, 
it  was  considerably  smaller  than  the  reading  accuracy  with  which  phase 
angle  could  be  determined  from  the  data,   and  thus  was  not  considered 
in  the  reduction  of  flight  test  data. 
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FLIGHT   TEST   PROCEDURE 

In  this  investigation,  frequency  response  was  measured  at  the 
single  flight  condition  of  5000  ft.   density  altitude,  129.6  mph  true  air- 
speed,  and  2870  lbs.   gross  weight.     Each  test  was  commenced  at 
steady  level  flight  with  the  amplidyne  control  set  at  zero  speed,   the 
oscillating  mechanism  exactly  centered  to  its  mid  position,   and  the 
lateral  controls  locked  in  this  trim  position.     The  run  was  then  initi- 
ated by  setting  the  speed  control  and  allowing  steady  state  conditions 
to  become  established.    The  necessity  for  driving  the  lateral  controls 
about  the  exact  trim  position  was  so  great  that  the  lock-on  system  to 
a  moving  oscillator,   as  employed  in  a  somewhat  similar  longitudinal 
study  (Ref.  1)  was  not  feasible.     In  spite  of  the  great  care  exercised 
to  center  the  oscillations  about  a  straight  and  level  flight  path,  prac- 
tically every  run  was  ultimately  terminated  by  the  excitation  of  the 
spiral  mode.     At  this  time  the  rapid  unlocking  features  in  the  mechan- 
ical system  design  were  most  appreciated. 

In  setting  the  speed  control,   an  approximation  to  the  desired 
frequency  could  be  obtained  by  the  use  of  a  voltmeter  across  the  ampli- 
dyne control  field  calibrated  to  read  frequency.     The  relatively  short 
duration  of  each  run  did  not  permit  fine  adjustments  to  obtain  a  specific 
frequency,   so  several  runs  around  a  desired  frequency  (primarily  near 
the  dutch  roll)  were  employed.     In  addition,   it  was  found  that  the   low 
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frequency  runs  were  best  started  by  going  rapidly  to  a  high  frequency, 
and  then  coming   down  to  the   desired  frequency.     Otherwise,   the 
response  to  the  slow  initial  deflection  drove  the  aircraft  so  far  into  the 
spiral  divergence  that  the  subsequent  control  reversal  was  insufficient 
to  establish  an  oscillation  about  a  straight  flight  path.     Although  it  had 
been  expected  that  all  controls  except  the  driven  one  would  be  locked 
during  a  run,  it  was  found  that  slight  elevator  movements  were  required 
to  maintain  a  constant  airspeed  throughout  the  run,  particularly  during 
oscillations  of  the  rudder . 
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RESULTS 

Data  Reduction 

The  quantities  measured  during  the  flight  test  program  were    &c\ 
or     <£r     ,     A     $     <p     »  and    1/s     .     These  four  quantities  were  teleme- 
tered to  the  ground  station,   recorded  on  magnetic  tape,   and  simulta- 
neously plotted  on  the  Sanborn  Recorder.     This  raw  data  had  noise  and 
distortion  superimposed  on  the  fundamental  sinusoidal  driving  frequency. 
An  example  of  such  a  record  is  shown  in  Fig.  11.     In  order  to  study  the 
response  of  the  aircraft  to  only  the  fundamental  of  the  driving  frequency, 
a  KROHN-HITE  Ultra-Low  Frequency  Band-Pass  Filter  was  used  to 
filter  both  the  driving  function  and  the  aircraft  responses.     From  the 
original  record  of  the  raw  data  the  fundamental  frequency  for  a  particu- 
lar run  could  be  determined  with  sufficient  accuracy  to  establish  the 
high  and  low  settings  for  the  band  pass  filter.     With  the  filter  thus    set 
to  pass  only  the  specific  frequency,   all  four  quantities  were  successively 
played  from  the  magnetic  tape  through  the  telemeter  receiver  (convert- 
ing the  pulse  width  coding  to  a  varying  voltage  output),  thence  through 
the  filter  to  the  Sanborn  recorder.     As  each  quantity  was  being  filtered 
and  recorded,   an    unfiltered  trace  of  the  forcing  function  was  also 
placed  on  this  new  record  to  act  as  a  time  reference.     The  results   of 
this  time  and  recorder  paper  consuming  process  were  four  sinusoidal 
records  of  the  measured  quantities,   each  with  a  common  unfiltered  time 
reference.     Fig.  12  shows  an  example  of  such  a  filtered  record  with  the 
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unfiltered  forcing  functions  matched  for  each  response.     From  these 
matched  responses  the  amplitude  ratios  and  phase  angles  relative  to 
the   driving  quantity  could  be   determined.      Since   the    settings   on  the 
filter  and  the  recorder  were  not  altered  during  the  filtering  of  a  given 
frequency,   any  amplitude  attentuation  or  phase  shifting  was   applied 
equally  to  all  four  quantities.     The  gain  of  the  recorder  was  initially 
set  for  each  frequency  to  accommodate  the  quantity  with  the  greatest 
amplitude,    such  that  this  quantity  would  utilize  almost  full  scale  on  a 
recorder  paper  channel.     By  using  this  technique,   and  increasing  the 
paper  speed  for  the  higher  frequencies,   a  rather  accurate  determina- 
tion of  the  amplitude  ratios  and  phase  angles  was  possible. 

The  quantities      Cj>    and      i^>  were  determined  from  the  records 
of      <p      and      ~)£     by  integration.,     This  was  readily  performed  by  divid- 
ing by     cO       and  subtracting  90°  from  the  phase  angle,    since  the   re- 
sponse was  sinusoidal.     The  magnitude  of  each  response  was  divided 
by  the  magnitude  of  the  forcing  function  and  the  appropriate  scale  fac- 
tors were  applied  to  give  the  desired  amplitude  ratios.     Since  the  sideslip 
vane  was  not  located  at  the  c.  g.   of  the  aircraft  an  additional  correction 
to  the     /$      response  was  required.     This  correction  was  a  function  of 
the  distance  of  the  sideslip  vane  from  the  e.g.  of  the  aircraft  and  is 
given  by  the  following  relation:     A  (8  =     —  .  OZ  7 4-    "¥~> 
The  final  frequency  responses  were  tabulated  in  Tables  II  and  in  and 
were  plotted  in  Figs.  13-18  and  labeled    "experimental." 
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Also  plotted  on  the  appropriate  graphs  in  Figs.  13-18  were  fre- 
quency response  curves  denoted  as    "theoretical."      These  curves  were 
the  result  of  an  initial  study  of  the  Navion  airplane  based  on  lateral 
stability  derivatives  determined  from  previous  tests,   such  as  described 
in  Ref.  2,   or  from  theoretically  determined  derivatives.     These  initial 
"best  guess"   derivatives  are  tabulated  in  Table  IV,  and  were  used  as  a 
rough  check  on  the  results,   as  well  as  the  basis  for  the  initial  system 
design. 

Determination  of  Derivatives 

Ref.   3  develops  a  method  for  extracting  stability  derivatives  from 
frequency  response  data  for  the  lateral  case.     The  method  involves  the 
application  of  the  vector  technique  to  the  three  linear  equations  of  lateral 
motion  to  separate  them  into  six  equations.     The  vector  responses  are 
resolved  parallel  and  perpendicular  to  the  control  deflection  vectors. 
The  perpendicular  resolution  eliminates  the  control  deflection  deriva- 
tives in  one  of  the  two  equations  derived  from  each  basic  equation  and 
thus  permits  solution  of  the  aerodynamic  derivatives  separately.     The 
development  of  these  equations  is  shown  in  Appendix  A.     The  method  of 
least  squares  is  then  applied  to  the  data  for  each  of  these  equations  to 
yield  lateral  stability  derivatives. 

Although  the  method  of  least  squares  is  probably  the  best  way  to 
eliminate  random  error,  it  is  not  readily  applicable  to  hand  calculation 
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techniques.     Automatic  computing  facilities  were  not  available  here, 
and  since  time  was  limited,   it  was  decided  to  use  the  simpler  method 
of  averages  in  conjunction  with  the  vector  equations  developed  in  Ref„  3 
to  determine  the  lateral  stability  derivatives  of  the   Navion.     The    K 
notation  of  Ref .   3  was  modified  to  a  starred  derivative  notation  for 
greater  clarity.     This  notation  is  defined  in  Appendix  A.     The  method 
of  averages,  as  described  in  detail  in  Ref.  4,  is  also  outlined  in 
Appendix  A. 

The  components  of  the  responses  along  and  perpendicular  to  the 
control  deflection  axis  are  defined  in  the  form  A+iB  with  the  appropri- 
ate subscripts.     Formulas  for  the  determination  of  the  responses  and 
their  derivatives  in  this  form  are  given  in  Appendix  A. 

The  two  equations  developed  from  the  basic  side  force  equation 
by  the  vector  technique  are  equations   (13)  and  (14)   in  Appendix  A. 
The   equation  describing  the  relation  between  the  components  perpen- 
dicular to  the  control  deflection  axis  was: 

(13)    Cy/g    -    CL    B<p    «     -«J   (A(2  +  Av) 

Since  the  only  unknown  was    Cy^*    it  was  possible,   in  theory,  to  deter- 
mine its  value  from  the  responses  to  either  an  aileron  or  rudder  forcing 
function  at  a  single  frequency.     Since  the  responses  in  yaw  and  sideslip 
were  greater  for  a  rudder  input,  it  was  felt  that  the  rudder  data  would 
give  a  better  solution  for    Cy^    •     In  order  to  eliminate  random  error 
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the  responses  at  12  selected  frequencies  were  taken  from  the  faired 
rudder  curves.     The  summation  of  the  12  equations  determined  from 
these  12  points  gave  one  average  equation,  from  which  the  value  of 
Cy^*    was  determined.     This  gave  a  value  of  Cy^  of  -0.591.     The  same 
technique  using  aileron  data  resulted  in  a  positive  value  for    Cyv,  , 
which  was,   of  course,  impossible.     Ref.   3  mentions  that  the  accurate 
determination  of    Cy^   may  be  difficult  from  this  equation,   since  it 
involves  the  sum  of  the  yaw  and  sideslip  components,   which  are  nearly 
equal  but  opposite  in  sign.     The  authors  of  Ref.   3  recommend  obtaining 
lateral  acceleration  data  and  using  a  modified  side  force  equation  to 
avoid  this  condition.     The  large  responses  in  yaw  and  sideslip  to  rudder 
inputs  tended  to  minimize  this  condition,  but  the  determined  value   of 
Cy^   was  nevertheless  subject  to  considerable  error. 

The  equation  describing  the  relation  between  the  components 
parallel  to  the  control  deflection  axis  and  therefore  including  the  deriva- 
tive,   Cy/    *  ,    was: 

(14)   Cy/A^  -  C*  A0   -   Cy<r*  +.  uj(se+  8-*) 

Using  the  previously  determined  value  for   Cy£     in  an  average  equation 
for  the  12  selected  points  yielded  a  value  for    Cy^  r    .      This  gave  a 
value  for    Cy^       of  0.206.     The  accuracy  of    Cy  ^       was  clearly  depend- 
ent on  the  correctness  of    Cy  .     ,   and  was,  in  addition,   subject  to  small 
differences  between  yaw  and  sideslip  components.     Cy  /     was  therefore 
of  doubtful  reliability. 
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The  vector  technique  gave  equations  (15)  and  (16)  of  Appendix  A 
from  the  basic  roll  equation.     The  relation  between  the  components 
perpendicular  to  the  control  deflection  axis  was: 

(15)  C«/b*  +  Cj>*ujA<p  4-  5g*>x&i,  -  C**u)Av  -  a>*54> 

Since   1^  was  very  small  compared  to  Ix    the  term  -g —     was  assumed 
to  be  zero.     Solution  for  the  three  unknowns  required  three  equations. 
The  12  frequency  points  from  the  aileron  data  were  now  grouped  to  in- 
clude four  frequencies  below  the  dutch  roll,  four  near  the  dutch  roll, 
and  four  above  the  dutch  roll,   and  average  equations  were  determined 
for  each  group  by  adding  the  appropriate  coefficients.     Solution  of  these 
equations  for    C].g*»     Cl    *,  and    Cj  *    were  converted  to    C|      , 
Ci    ,    and    Ci    ,     which  were  -0.0505,    -0.436,  and  0.0390,  respectively. 
In  the  solution  of  these  equations  it  was  found  that  an  accurate  determin- 
ation of    Ci.    was  difficult,    since  the  solution  involved  a  small  differ - 
r 

ence  between  large  quantities;   and  small  errors  were  magnified  many 
times.     An  attempt  was  made  to  carry  greater  accuracy  in  the  solution 
than  the  data  actually  warranted,   but  of  course  this  procedure,   while 
required  to  give  any  solution  at  all,   could  not  give  an  accurate  solution 
for    Cir.      This  condition  was  not  surprising,   since  it  is  well  known 
that  the  cross  derivatives  are  difficult  to  determine  from  this  sort  of 
data. 

Solution  of  the  same  equation  using  a  grouping  of  points  from  the 
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rudder  data  should  have  given  the  same  derivatives  as  previously 
determined  from  the  aileron  data-     When  this  solution  was  attempted, 
however,   it  was  found  that  the  simultaneous  equations  were  ill-condi- 
tioned.    This  arose  from  the  fact  that    (3     and    <p    varied  in  nearly 
the  same  manner  over  the  entire  frequency  range,   both  in  phase  angle 
and  magnitude.     Additional  frequency  points  were  taken  and  various 
groupings  were  tried  in  an  attempt  to  eliminate  this  difficulty,  but  all 
were  unsuccessful.     Ref.   3  suggests  a  possible  method  of  avoiding  the 
difficulty  by  eliminating  the  cross  derivative  coefficient  by  resolution 
of  the  vectors  along  the    (fj     axis,   instead  of  along  and  perpendicular 
to  the  control  deflection  axis.     This  was  attempted,  but  solution  of  the 
modified  equation  was  still  impossible. 

B:  was  felt  that  the  rudder  data  was  slightly  better  than  the  aileron 
data,  due  to  refinements  in  the  flight  test  technique  as  the  program  pro- 
gressed,  and  it  was  therefore  highly  desirable  to  include  the  rudder 
data  in  the  solution  for  the  roll  derivatives.     The  solution  of  the  aileron 
equations  gave  a  value  of    Ci         which  was  felt  to  be  quite  accurate, 
since  it  was  in  close  agreement  with  theory  and  previous  experimental 
results.     This  value  was  assumed  known  in  the  rudder  equations  for  the 
rolling  motion,  and  the  two  remaining  derivatives    Ci6       and    Ci 
were  determined  by  taking  two  groups  of  six  frequency  points  each  and 
arriving  at  two  equations  in  the  two  desired  unknowns .     The  values  de- 
termined for    C\        and    C\       were  -0. 0513  and  0.119  respectively. 
p  * 
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It  was  interesting  to  note  that  the  value  of    C\g    determined  in  this 
manner  was  nearly  the  same  as  that  determined  from  the  aileron  data, 
and  the  doubtful  derivative,    Cir,    was  nearer  the  theoretically  deter- 
mined value.     The  derivatives  determined  in  this  manner  were  arbi- 
trarily chosen  as  the  final  values,    since  the  authors  felt  that  they  repre- 
sented the  best  possible  combination  of  all  the  experimental  data. 
The  relation  between  the  components  parallel  to  the  control 
deflection  axis  was: 

(16)  a*   A*-   C4u>z<p  +    ^    ^Af  +    C^wB^-C**  -  ^A<J> 


.fi     r^(i  v.<^    cu  ocp    -r         j 


< 


Using  the  previously  determined  values  of    Cj    *,     Cj   *,    and    Ci   * 
and  utilizing  one  average  equation  based  on  all    12  frequency  points 
from  the  rudder  data  permitted  solution  for    C\,     *,     which  gave  a 
Cj  r        of  0.  0128.     Similarly,   solution  of  the  same  equation  using  the 
aileron  data  gave    Cl^     *,    which  resulted  in  a  value  for    Ci^      of 
-0.108. 

Expressing  the  basic  yaw  equation  in  vector  form  resulted  in 
equations   (17)   and  (18)   of  Appendix  A.     For  the  components  perpen- 
dicular to  the  control  deflection  axis  the  relation  was: 


(17)    -Cn*  %   -  C^   *o4<*  -     v*  ^^  +cnN  A^=-    <oaB> 


Because  of  the  basic  similarity  between  the  equations  for  the  yawing 
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and  rolling  motions,   the  unknowns  in  equations    (15)   and  (17)   have  the 
same  numerical  coefficients,   with  only  the  constant  term  differing. 
Because  of  this,   rapid  solution  of  the  three  average  equations  for  the 
yawing  motion  was  possible  using  the  same  grouping  of  frequency 
points  from  the  responses  to  aileron  oscillations.     This  gave  values  of 
CnA    »     C7i      »    and    cn.r      °*  0.0771,    -0.129,   and  -0.0841,   respectively, 
The  same  difficulty  existed  here  in  the  determination  of    Cnp    as  had 
previously  existed  in  the  determination  of    C\r. 

Because  the  numerical  values  of  the  coefficients  of  equation  (17) 
were  the  same  as  those  of  equation  (15),  the  three  simultaneous  equa- 
tions for  the  yawing  motion  derived  from  the  rudder  data  were  also  ill- 
conditioned  for  solution.     The  method  applied  to  solve  equation  (15) 
using  an  apparently  accurate  value  of    Cj      could  not  be  applied  in  this 
case,  because  both    C^^      and    Cn^      as  determined  from  the  aileron 
data  were  felt  to  be  doubtful,   since  they  were  not  in  agreement  with 
either  theory  or  the  results  contained  in  Ref.   5.     In  this  case  final 
solution  was  effected  by  solution  of  a  set  of  simultaneous  equations, 
two  of  which  were  determined  from  the  aileron  data  at  frequencies  near 
and  above  the  dutch  roll,   and  the  third  representing  the  average  of  all 
the  rudder  data.     This  approach  appeared  legitimate  since  each  equa- 
tion in  the  set  was  independent  of  the  type  of  forcing  function.     This 
particular  choice  of  equations  resulted  from  the  opinion  of  the  authors 
that  the  aileron  data  was  better  at  higher  frequencies  and  that  inclusion 
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of  the  generally  more  reliable  rudder  data  was  desirable.     It  should  be 
emphasized  that  this  particular  choice  was  peculiar  to  this  experiment. 
Solution  of  these  equations  gave  values  for    Cn^    ,     C?ip    •   an<^  ^nr 
of  0.0895,    -0.106,   and  -0.0908  respectively. 

The  equation  relating  the  components  parallel  to  the  control 
deflection  axis  for  the  yawing  motion  was: 

(18)-Cn/A/S  +  Cn*  lo&$+XJ£  «,M<t>  -  Cn*  lu^.  -  C^  =  u^A* 

Using  the  previously  determined  values  of    Cn      ,     C^   »    and    Cnr 
and  utilizing  all  12  frequency  points  from  the  rudder  data  permitted 
solution  for    O^  e     *,    which  gave  a    C^  ,      of   -0„0673.     Similarly  solu- 
tion of  the  same  equation  using  the  aileron  data  gave    Cn  tA  *»    which 
resulted  in  a  value  for    Cn  <       of  0.00376. 

The  final  derivatives  determined  in  this  investigation  are  tabulated 
in  Table  IV.     This  table  also  contains  the  original  theoretically  deter- 
mined derivatives,  and  the  derivatives  determined  for  the  same  aircraft 
by  an  analog  computer  solution  as  outlined  in  Ref.   5. 

In  order  to  check  the  accuracy  of  the  method  used  to  determine  the 
aircraft's  stability  derivatives,  the  transfer  function  coefficients  based 
on  them  were  determined  as  outlined  in  Appendix  A,   and  from  these  the 
frequency  responses  were  calculated.     The  experimental  transfer  func  - 
tions  are  tabulated  in  Table  V,  and  the  frequency  responses  based  on 
them  are  shown  in  Fig.  19-24. 
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DISCUSSION 

Dynamic  flight  testing  is  becoming  increasingly  important  with  the 
advent  of  high  speed  high  altitude  aircraft  because  the  aircraft  designer 
is  more  concerned  than  ever  before  with  the  accurate  determination  of 
stability  derivatives  and  the  associated  transfer  functions.     Although 
transient  methods  are  the  most  popular  at  the  present  time,   it  was  felt 
that  if  the  instrumentation  requirements  of  steady  state  dynamic  flight 
testing  could  be  reduced  and  the  necessity  for  absolutely  sinusoidal  con- 
trol inputs  without  the  added  requirement  of  harmonic  analysis  could  be 
eliminated  the  forced  oscillation  technique  for  determining  frequency 
response  and  thus  stability  derivatives  and  transfer  function  coefficients 
could  still  be  used  to  advantage  under  certain  circumstances.     It  was  the 
purpose  of  this  investigation  to  examine  the  possibility  of  obtaining  ade- 
quate frequency  response  with  relatively  simple  instrumentation,   and  if 
possible  to  eliminate  the  necessity  of  harmonic  analysis  due  to  distorted 
sinusoidal  inputs. 

Since  investigation  of  the  methods  involved  was    of  primary  con- 
cern here,  flight  tests  were  conducted  at  only  one  altitude,  airspeed, 
and  gross  weight  condition,   although  it  should  be  emphasized  that  in  any 
actual  flight  test  program  to  determine  derivatives,   all  of  these  condi- 
tions should  be  varied  over  the  operating  range  of  the  aircraft. 

Before  the  design  of  the  instrumentation  system  could  be  accom- 
plished a  prediction  of  the  results  based  on  theoretically  determined 
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derivatives  and  control  gradients  was  required.     Obviously  the  experi- 
mental data  had  to  include  frequencies  on  either  side  of  the  dutch  roll 
with  frequencies  at  the  dutch  roll  being  of  primary  interest.     Previous 
investigations  of  the  lateral  dynamic  stability  of  the  Navion  indicated 
that  the  period  of  the  dutch  roll  was  between  2   and   3   seconds.     An 
analysis  based  on  the  best  available  derivatives,   some  calculated  theo- 
retically from  Ref.   6  and  others  taken  from  Ref.  2  and  7  was  performed. 
This  indicated  that  the  period  of  the  dutch  roll  was  about  2.4  seconds. 
Theoretical  frequency  responses  to  both  rudder  and  aileron  forcing 
functions  were  calculated  and  these  served  to  design  the  instrumentation 
system. 

It  was  decided  to  design  the  oscillator  to  cover  frequencies  be- 
tween 0.  5  and  10  radians  per  second.     The  upper  limit  proved  higher 
than  necessary  and  better  results  would  have  been  obtained  if  the   design 
had  limited  the  upper  frequency  to   5  radians  per  second.     This  upper 
limit  was  found  to  be  well  above  the  dutch  roll  frequency  and  oscilla- 
tions at  higher  frequencies  produced  only  negligible  responses .     Exam- 
ination of  the  control  deflections  necessary  to  produce  recordable 
responses  without  exceeding  the  linear  limit  of  the  equations  and  the 
control  gradients  led  to  the  decision  to  use  a    ^  H.P.  motor  and  an 
amplidyne  to  control  its  armature  current.     Early  experimentation  with 
a  mechanical  variable  reduction  gear  indicated  that  this  method  was 
unsatisfactory. 
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Ref .  1  recommends  a  careful  preliminary  study  of  the  dynamics 
of  the  recording  instruments  since,   if  their  natural  frequency  is  not 
much  larger  than  the  natural  frequency  of  the  aircraft,  large  errors 
arise  which  must  be  corrected  for.     This  was  done,  and  it  was  found 
that  the  natural  frequency  of  the  rate  gyros  and  sideslip  vane  were  such 
that  the  attenuation  of  their  signal  outputs  was  effectively  zero  and  that 
the  phase  shift  introduced  was  less  than  0.5°.     Since  other  errors  in 
phase  angle  measurement  were  expected  to  be  of  considerably  greater 
magnitude,  no  attempt  was  made  to  correct  for  the  dynamic  response 
of  the  recording  instruments. 

The  instrumentation  system  based  on  the  results  of  the  above  cal- 
culations proved  completely  adequate  with  the  exception  of  the  oscillator. 
It  was  found  that  since  the  output  torque  was  directly  proportional  to  the 
armature  current  of  the  motor  the  torque  at  low  frequencies  and  thus 
low  armature  currents  was  insufficient  to  produce  anything  approaching 
a  pure  sinusoidal  control  motion.     Consideration  was  given  to  changing 
the  reduction  gear  ratio  and  thus  providing  more  torque  at  low  frequen- 
cies, but  this  was  discarded  because  of  time  requirements.     The  gyros 
were  properly  utilized,   and  their  limitations  were  never  exceeded;   and 
the  original  estimates  of  the  control  deflections  required  to  produce  the 
desired  responses  proved  accurate,  making  major  modification  of  the 
original  instrumentation  unnecessary. 

Several  difficulties  were  encountered  during  the  flight  test  program. 
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The  initial  flights  were  conducted  with  the  rudder  locked  and  with  the 
ailerons  driven  by  the  oscillator.,     At  oscillation  frequencies  well 
above  the  dutch  roll  the  steady  state  condition  was  rapidly  reached  and 
satisfactory  data  could  be  obtained  before  the  spiral  mode  was  excited. 
At  frequencies  near  the  dutch  roll  the  spiral  mode  was  frequently  ex- 
cited before  the  steady  state  condition  was  reached  necessitating  ter- 
mination of  the  run  before  any  data  could  be  obtained.     At  frequencies 
below  the  dutch  roll  the  initial  control  deflection  was  sufficient  to  ex- 
cite the  spiral  mode^and  no  data  at  all  could  be  obtained.     After  some 
experimentation  it  was  found  that  most  of  this  difficulty  could  be  elim- 
inated by  starting  the  oscillations  at  a  high  frequency  and  then  reducing 
the  frequency  to  the  desired  value,   and  this  procedure  was  adopted  in 
all  subsequent  flight  tests  both  for  aileron  and  rudder  oscillations. 
During  the  runs  involving  aileron  oscillations  there  was  very  little  indi- 
cation to  the  pilots  when  the  dutch  roll  frequency  was  being  approached, 
and  it  was  not  until  the  data  was  reduced  that  it  could  be  seen  that  the 
responses  peaked  at  this  frequency.     During  the  runs  involving  rudder 
oscillations,   however,  the  responses  became  quite  violent  at  the  dutch 
roll  and  apparently  very  sensitive  to  airspeed.     During  the  first  rudder 
run  at  the  dutch  roll  frequency  the  airspeed  control  was  poor,   and  as 
speed  decreased  by  as  little  as   6  mph  the  responses  decreased  in  a 
very  noticeable  manner.     As  speed  was  increased  to  the  desired  value 
by  retrimming  the  elevator  ^  the  resonant  condition  was  again  reached 
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and  the  responses  became  large.     A  similar  variation  was  noted  with,  an 
increase  in  airspeed  over  the  desired  value,   but  the  effect  was   not   as 
pronounced.     The  difference  between  the  effects  was  explained  by  refer- 
ence to  the  final  response  curves  which  show  a  rather  flat  top  at  the 
dutch  roll  frequency.     A  decrease  in  velocity  corresponds  to  an  increase 
in    T'    and  this  in  turn  corresponds  to  a  decrease  in  the  dutch  roll  fre- 
quency.    Since  the  actual  frequency  of  the  oscillations  was  maintained 
constant,  the  operating  point  at  the  decreased  velocity  was  actually 
above  the  dutch  roll.     Since  the  experimentally  determined  response 
curves  fall  off  more  rapidly  above  the  dutch  roll  than  below  it,   the  con- 
dition was  satisfactorily  explained.     Because  of  the  violence   of  the 
responses  near  the  dutch  roll  and  their  sensitivity  to  airspeed  great 
care  was  necessary  to  keep  the  speed  at  its  desired  value  and  runs 
where  speed  control  was  poor  were  discarded. 

Ref.  1  discusses  some  difficulty  in  the  design  of  satisfactory  con- 
trol locking  and  trim  devices  for  the  oscillation  mechanism.     The  orig- 
inal design  accomplished  here  utilizing  the  recommendations  of  Ref.  1 
proved  highly  successful  and  no  difficulty  was  encountered  either   in 
obtaining  proper  trim  or  in  releasing  the  oscillator  or  the  locked  control 
surface  when  the  attitude  of  the  aircraft  made  this  necessary. 

While  it  was  felt  in  the  original  design  that  a  measure  of  the  con- 
trol field  voltage  to  the  amplidyne  would  serve  as  a  satisfactory  indica- 
tion of  the  frequency  of  oscillation,  it  was  found  that  this  assumption 
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was  not  wholly  justified.     A  rough  determination  of  frequency  could  be 
made  from  this  voltage  reading,  but  because  of  residual  magnetism  in 
the  control  field  of  the  amplidyne,   the  desired  accuracy  could  not  be 
attained.     Precise  frequency  determination  for  data  analysis  was  made 
from  the  recorded  traces 5but  since  a  fairly  accurate  knowledge  of  fre- 
quency in  the  cockpit  during  the  flight  testing  was  desirable  it  is  recom- 
mended that  in  subsequent  tests  of  this  nature  a  tachometer  generator 
be  incorporated  on  the  driving  motor. 

In  the  analysis  of  the  data  at  low  frequencies  it  was  found  that, 
although  the  inputs  were  not  purely  sinusoidal,  the  steady  state  condition 
was  rapidly  reached.     Harmonic  analysis  of  both  the  responses  and  the 
forcing  function  was  considered  but  was  discarded  as  being  too  time  con- 
suming.    Since  all  the  data  had  been  recorded  on  magnetic  tape,  a  tech- 
nique of  filtering  both  the  forcing  function  and  the  aircraft  responses 
with  a  low  frequency  band-pass  filter  was  developed.     The  procedure,  as 
previously  described,  was  really  nothing  more  than  a  fast  and  simple 
method  of  performing  harmonic  analysis  of  the  traces  and  required  that 
the  data  be  recorded  on  magnetic  tape.     The  technique  was  a  valuable 
one,  and  it  was  felt  by  the  authors  that  it  could  be  utilized   to  advantage 
in  all  forms  of  steady  state  dynamic  flight  testing  where  telemetering 
and  magnetic  recording  of  data  is  performed. 

The  use  of  telemetering  equipment  in  a  flight  test  program  of  this 
sort  requires  more  complicated  instrumentation  than  the  simple  SFIM 
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recorder  utilized  in  Ref .  1,   in  that  a  ground  station  and  its  associated 
equipment  is  required.     The  actual  instrumentation  in  the  aircraft  is 
little  complicated  by  the  substitution  of  the  telemeter  transmitter   for 
the  in-flight  recorder.     The  utility  derived  from  having  magnetic   tape 
recordings  when  the  telemeter  scheme  is  used  is  enormous,   and  it  is 
felt  that  the  advantages  of  this  method  outweigh  the  complications  pre- 
sented by  the  extra  instrumentation  required,,     The  filter  technique   of 
course  would  be  impossible  without  tape-recorded  data,  but  in  addition 
tape  records  of  this  sort  may  be  expanded  to  any  desired  time    scale 
making  phase  angle  measurements  easier  and  more  accurate.     It  is 
also  possible  to  amplify  small  responses  to  distinguishable  values,  but 
care  must  be  taken  not  to  exceed  the  accuracy  of  the  record.     The  tele- 
metering scheme  utilized  here  was  accurate  to  1%  of  full  scale  and  thus 
small  responses  could  be  expected  to  have  a  much  greater  percentage 
error  than  large  ones.     It  was  therefore  desirable  during  this  investiga- 
tion to  obtain  responses  which  were  as  large  as  possible  without  exceed- 
ing the  range  of  linearity  of  the  equations  of  motion.     This  was  accom- 
plished by  varying  the  forcing  function  amplitudes  over  the  frequency 
range  to  give  desirable  responses.     In  general,   large  control  deflections 
were  required  at  high  frequencies  and  small  ones  at  and  below  the  dutch 
roll  frequency. 

The  method  used  in  this  investigation  to  determine  stability  deriva- 
tives from  frequency  response  was  basically  sound „     As  mentioned 
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previously,    considerable  difficulty  was  encountered  in  the  determination 
of  the  cross  derivatives,   and  it  was  felt  that  accurate  determination  of 
these  would  require  much  more  accurate  frequency  response  data  than 
could  be  obtained  with  the  simple  instrumentation  utilized  here.      The 
accurate  determination  of  phase  angles  from  data  obtained  using  the 
forced  sinusoidal  oscillation  technique  is  difficult  indeed,   and  phase  angle 
errors  can  have  a  large  effect  in  the  equations  for  the  solution  of  the 
stability  derivatives,   especially  when  the  coefficients  of  these  equations 
involve  sines  of  angles  near   0°  or   180°  or  cosines  of  angles  near   90° 
or   270°o     The  technique  developed  here,   which  involved  fairing    a 
smooth  curve  through  the  experimental  points  and  taking  data  for  the 
computation  of  the  derivatives  from  the  smooth  curve,   eliminated  some 
of  this  difficulty,   although  it  was  felt  that  the  method  of  averages  could 
probably  eliminate  most  of  this  random  error  even  if  individual    data 
points  were  used. 

The  equations  derived  from  the  rudder  data  were,   as  previously 
mentioned,   ill-conditioned  for  solution,   and  this  fact  required  determin- 
ation of  some  of  the  derivatives  from  the  aileron  data  before  the  rudder 
equations  could  be  used.     In  theory,   it  is  possible  to  determine  all  the 
derivatives,   except  of  course  the  aileron  control  derivatives,   from 
rudder  data  alone „     While  this  was  not  possible  for  the  Navion  aircraft 
at  the  condition  investigated,   it  was  felt  that  for  other  aircraft  or  other 
flight  conditions  the  difficulties  encountered  here  might  not  be  present. 
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It  became  apparent  to  the  authors  that  a  need  exists  for  extensive  study 
of  the  determination  lateral  stability  derivatives  from  frequency  re- 
sponse data,   especially  as  concerns  ill-conditioned  equations  and  possi- 
ble methods  for  their  solution.     While  it  finally  proved  possible  to 
obtain  fairly  accurate  derivatives  in  this  investigation,   certain  devia- 
tions from  the  basic  method  were  required.     The  methods  applied  here 
were  peculiar  to  this  particular  aircraft  and  flight  condition  and  might 
very  well  prove  unsuccessful  for  other  aircraft  or  flight  conditions.     The 
best  procedure  to  follow  in  the  extraction  of  stability  derivatives  from 
frequency  response  data  is,   at  the  present  time,   a  function  of  the  type 
and  accuracy  of  the  data,   although  the  authors  feel  that  with  time  a 
standardized  procedure  for  eliminating  difficulties  could  be  developed. 

The  use  of  the  method  of  averages  for  eliminating  random  error 
instead  of  the  method  of  least  squares,   as  outlined  in  Ref.  4,  proved 
quite  satisfactory  and  had  the  added  advantage  of  reducing  the  computa- 
tion time  considerably.     When  using  this  method  care  must  be  taken 
when  grouping  the  points  for  averaging.     If  automatic  computing  facili- 
ties are  available  and  large  amounts  of  data  require  reduction,    the 
method  of  least  squares  would  probably  prove  more  satisfactory  and 
less  subject  to  error. 

It  was  difficult  in  this  analysis  to  acquire  a  feel  for  the  accuracy 
of  many  of  the  individual  derivatives.     As  previously  stated,  the  nature 
of  the  equations  and  the  fact  that  small  differences  of  large  quantities  > 
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were  involved  in  the  solution  for    Cy     ,     Cir»    and    CL      made  the 
accuracy  of  these  derivatives  questionable.     In  the  solution  for    Cy        , 
C\,      ,    and    Cn  ,        the  controlling  terms  involved    Cy^  ,     Ci   ,     and 
C™      respectively,   and  thus  these  control  derivatives  were   also   of 
doubtful  accuracy.     Lateral  acceleration  data  would  have  eliminated 
this  difficulty  in  the  case  of  the  side  force  derivatives,  but   was   not 
obtained  here.     In  the  determination  of  the  other  derivatives,  however, 
the  potential  degree  of  accuracy  was  not  evident  from  the  equations. 

A  solution  of  the  characteristic  equation  using  the  experimentally 
derived  derivatives  resulted  in  a  period  of  2.  38  seconds  and  a  time  to 
half  amplitude  of  1.69  seconds  for  the  dutch  roll  mode.     Transient  data 
obtained  for  the  same  aircraft  at  the  same  flight  conditions  indicated  a 
period  of  2.40  seconds  and  a  time  to  half  amplitude  of  1.70  seconds. 
In  view  of  this  close  agreement,   it  was  concluded  that    Cn^  .    which 
primarily  determines  the  period  of  the  dutch  roll  mode,   and    Cn   ,  the 
primary  damping  derivative,   were  accurately  determined. 

A  comparison  of  the  frequency  responses  derived  from  the  experi- 
mentally determined  derivatives  and  the  experimental  points  indicated 
certain  regions  of  discrepancy.     In  Fig.  19  the  calculated  amplitude 
response  in  roll  to  aileron  oscillations  was  noticeably  low  at  and  above 
the  dutch  roll  frequency,   and  the  calculated  phase  angle  was  slightly 
high  in  the  same  region.     An  analysis  of  the  transfer  function  contained 
in  Table  V  indicated  that  an  increase  in  the  coefficient  of  the  first  power 
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term  in  the  numerator  would  remedy  this  condition.     The  most  obvious 
means  of  accomplishing  this  was  an  increase  in  the  value  of    C^, 
An  increase  in  this  derivative  would  also  increase  the  constant  term  in 
the  numerator  of  the  yaw  transfer  function  and  thus  increase  the  ampli- 
tude response  at  low  frequencies.     Fig*   20  shows  the  need  for  this  cor- 
rection.    While  a  detailed  analysis  of  the  effect  of  increasing    C]^ 
was  not  performed,  these  considerations  led  to  the  conclusion  that  the 
value  determined  for    Clj       was  low,   and  a  comparison  with  the  results 
of  Ref.   5  indicated  this  to  be  the  case. 

In  Fig.  22  the  calculated  amplitude  response  in  roll  to  rudder 
oscillations  is  substantially  less  than  the  experimental  response  at  low 
frequencies  and  the  phase  angle  is  consistently  larger  for  all  frequencies 
An  analysis  of  the  transfer  function  for  this  response  indicated  that  an 
increase  in    Cn  ,        would  tend  to  correct  both  these  discrepancies.     In 
addition,  the  calculated  low  frequency  amplitude  response  and  high  fre- 
quency phase  response  for  sideslip,   as  shown  in  Fig„  24,   would  also  be 
altered  favorably  by  a  larger    Cn<£      .     A  comparison  with  the  results 
of  Ref.   5  indicated  that  the  value  determined  here  for    Cn  ,        was  too 
low. 

Since  the  match  between  the  calculated  and  actual  responses  was 
generally  good  with  the  exception  of  the  cases  mentioned  above,  it  was 
felt  that  the  determination  of  those  derivatives  with  the  greatest  effect 

on  the  frequency  response,   i.e.,     Ci      ,     Cn      ,     Ci    ,     and    C«    ,    was 

x  $  ■"#  p  r 
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good.     While  detailed  curve  fitting  of  individual  responses  by  altering 
individual  derivatives  was  not  considered  within  the  scope  of  this  inves 
tigation,    such  a  technique  for  correcting  derivatives  might  warrant 
further  study. 
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CONCLUSIONS  AND  RECOMMENDATIONS 

The  method  utilized  in  this  investigation  to  determine  the  lateral 
frequency  response  of  an  airplane  proved  entirely  feasible  and  could 
undoubtedly  be  applied  to  many  current  production  aircraft.  The 
method,  however,  had  the  disadvantage  of  requiring  more  flight  time 
than  currently  popular  transient  response  methods,  and  thus  would  not 
be  practical  for  missile  testing., 

Accurate  determination  of  frequency  response  previously  required 
an  elaborate  system  for  the  generation  of  an  exact  sinusoidal  forcing 
function  or  time  consuming  harmonic  analysis  of  the  datao     The  filtering 
technique  developed  here  eliminated  these  difficulties,   and  could  have 
wide  application  in  this  type  of  flight  testing. 

The  method  utilized  here  for  the  determination  of  stability  deriva- 
tives from  frequency  response  data  proved  accurate  for  those  derivatives 
with  significant  effect  on  the  frequency  response.     It  was,  however,   im- 
possible to  determine    On    »     Ci     ,     Cn  /      ,    and    Cl  ,         with  a  high 
degree  of  accuracy.     Better  results  could  have  been  realized  in  the 
determination  of    Cy       and    Cy  ,         if  lateral  acceleration  data  had  been 
obtained,   and  it  is  recommended  that  such  measurements  be  made  in 
future  tests.     It  is  further  recommended  that  a  study  be  made  to  develop 
a  systematic  method  for  eliminating  the  difficulties  encountered  here 
with  ill-conditioned  equations. 
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TABLE  I 


PHYSICAL  CHARACTERISTICS   CF  NAVION 


A.         WINQ 

1.  Total  Wing  Area,  S*  (includes  flaps,  ailerons 
and  19 087  Ft?"  covered  by  the  fuselage) 

2o  Span  by 

3.  MAC 

U*    Angle  of  Incidence 
Root,   ir 

Tip,  it 

5.  Twist 

Aerodynamic 
Geometric 

6.  Airfoil  Section 

Root 
Tip 

7o  Aspect  Ratio,  AR^ 

80  Taper  Ratio,     A  v 

9.  Dihedral 

10*  Root  Chord 

llo  Tip  Chord 


l8lio337  Ft.2 
33o378  Ft. 
68,352  in. 


♦2° 
-10 


2°  31 » 
3° 


NACA  Uil^R 
NACA  6U0R 

6>6kk 

0.5265 

7o5° 

7o2    Ft. 

3o92  Ft. 


A-lo     AILERON  (one  Aileron  only) 

(Beaded  skin,  Frise  Nose  Balance, 
No  Trailing  Edge  Extension) 

1.  Area,      Sa 

20  Span,     b& 

3*  Aileron  Deflection,    ^  a 

bo  Boost 

5.  Trim  Tab 


120O-120O 

Wheel  Throw 
Fixed  bend  tab  on 
right  aileron. 

20161Ft.2 

61.987  in. 

30°  UP,  20°  DN. 

None 

Fixed  bend  tab  on 
right  aileron  o 


B.         HORIZONTAL  TAIL 

1.  Total  Area 

(Includes  2,368  Ft,2  covered  by  fuselage) 

2,  Span,   bjj 

3»  Airfoil  Section 

U,  Taper  Ratio,    A  t 

5o  MAC,       Ct 

6,  Aspect  Ratio,      /5^  ^ 

B-lo     HORIZONTAL  STABILIZER 

1.  Area,  Se 

2,  ^g  (with  reference  to  FRP) 

B-2.     ELEVATORS 

(No  Trailing  Edge  Extensions) 


I,  Total  Area,     Se 

2  o  Span,     be 

3,  Deflection,     ^  e 

Uo  Trim  Tabs  (32»  Span,  U-l/2"  Chord) 

5o  Root  Chard 

6.  Tip  Chord 


13.051  Ft.2 

13ol72  Ft, 
NACA  0012 
.67 

3.3U  Ft. 

In  02 

28.953  Ft.2 

Smooth  Skin 
Flat  sided 
No  triai  bungee 
Balance  spring 

Hw098  Ft.2 
73.582  in. 
30°  UP,  20°  DN. 
30°  UP,  30°  DN. 

1.5  Ft. 

1.0  Ft. 


C.         VERTICAL  TAIL 


1.     Total  Area,     S, 


(Includes  2o577  Ft.2  blanketed  by 
fuselage  and  excluding  LI483  Ft.2 
of  dorsal  fin) 


Airfoil  Section 
Root 
Tip 


12.925  Ft,2 


NACA  0013 o2  Mod. 
NACA  0012 -6U  " 


C-l.     VERTICAL  STABILIZER 
1*     Area,     Sy 

6 

(Includes  0.1U27  Ft. 2  blanketed  by  fuselage) 
2. -X  (With  reference  to  fuselage  0) 

C-2.     RUDDER 


1.  Area,     Sj. 

2*  Rudder  Deflection,    ^  r 

3o  Rudder  Pedal  Throw 

Uo  Trim  Tab 


6.873  Ft.2 

2°  Nose  left 

Smooth  Skin 
Rigged  3°  Rt.  to 
Fin  0,    Fixed 
Bend""Tab. 

6.052  Ft.2 
17°  L  -  23°  R 

5.75  in. 
Fixed  Bend  Tab 


D.         PROPELLER 

lo  Activity  Factor,  A.F 

2.  Diameter,  d 

3.  Pitch  G  @  .75  R 


Hart »e  11 
100 
86" 
21©  30' 


Eo    NOSE  OSAR 
1*  Travel 
2o  Strut  Angle  (To  Vertical) 


♦  30° 
10° 


F. 


MISCELLANEOUS 


lo     Length,   overall 
20     Tail  length,     A  VT 


27o25  Ft. 
16.88  Ft. 


0. 


POWER 


1.     Engine 


Continental  E-185  Engine 
185  HP  @  2300  RPM,  29"  Hg 
M.P.  @  Sea  Level 


TABLE  H 

EXPERIMENTAL  FREQUENCY  RESPONSE  DATA 
TO  SINUSOIDAL  AILERON  OSCILLATIONS 


T.A.S. 

Density  Alt< 
Gross  Wto 


129o6  m.p.h. 
5000     ft. 
2870     lbs. 


60 

id) 

*♦*. 

1  <W 

1 

if  I 

Oo73 

3»20 

1 
92.0 

1.225 

i 

-  11.0 

.362 

173  o3 

loll* 

2d75 

76.0 

.703 

-    7o5 

1 

<»li85 

15202 

1.61 

1.1452 

88o2 

.637 

-  28o9 

o566 

liiO.O 

1.78 

1.322 

B6.0 

.601 

-  30.0 

.529 

Hi5oO 

2.23 

1.01*7 

*.i 

.8W1 

-51.0 

.920 

1080I 

2.33 

1.01*2 

1 

&.0 

.881 

-58.0 

08O6 

119.0 

2.57 

1.197 

8U.0 

i 

.790 

-105.0 

1.032 

72.8 

2o59 

1.U58 

88.5 

• 

.750 

-132.2 

1.005 

16. 0 

2.62 

lol30 

loli.U 

O820 

-  9U.0 

lolOO 

85  oO 

2.79 

1.UU8 

90.0 

.725 

-136.7 

.982 

U60O 

- 

2.97 

1.338 

83.2 

.666 

-132 08 

.855 

28.9 

: 

3.27 

1.001 

69.li 

o259 

+180.0 

.367 

8.7 

3.79 

J 

.899 

66.1 

.213 

172.0 

.283 

10.0 

3o92 

J         .822 

I       65.3 

0I6O8 

166  06 

.260 

356.5 

Uo63 

.681 

1 

63  o5 

.0851 

151.0 

■ 

5.20 

.562 

52.0 

.0616 

131  o3 

.1265 

33U.7 

5.B5 

.196 

59.0 

.OI1O8 

125.0 

.0775 

336o3 

T.A.S. 
Density  Alt, 
Gross  Wto 


TABLE  III 

EXPERIMENTAL  FREQUENCY  RESPONSE  DATA 

TO  SINUSOIDAL  RUDDER  OSCILLATIONS 


129.6  m.poh. 
5000     ft. 
2870     lbso 


1.05 

1.315 
j   1.615 
1.85 

1.87 

■ 

2o21 
2ol*l 
2.1*8 
2.58 
2o58 
2.69 
2.99 
3.71 

UoOO 

U.35 
l*o68 
U.70 
U.Q3 

i    5o3U 

I    5«38 

5.71 


7.21 


1.222 

1.121 

1.100 

1.106 

1.352 

lo358 

loU65 

1.U8 

1.55 

1.51 

1.5147 

loOO 

.152 

.291* 
.253 
.181 
.161* 
.160 
.101* 

.083 
.01*1*7 


100  o2 

914.5 
9i*.6 
79.1* 
82.5 
57.1 
37o3 
5Uo6 
21o6 

35ol4 

♦  5.2 
-2lu5 

-51.8 
-63.6 

-58.7 
-66.6 

-66  o5 
-68.0 
-77.6 

-75o5 
-83.5 


Sr 

7.U 

.901* 

1.053 
lo211 

1.3H 

I08IO 

2.165 

2ol3 

2.21 

20168 

20218 

I066I 

.822 

.599 

.1*93 

.390 

.36h 

.338 

.253 

.21*35 

.216 

.127 


176.2 
168  ok 
171.8 
160.0 
161*. 0 

U*1.0 

.119.9 

132  o5 

106  0I4 

12i*.o 

79o5 

73.2 

38.9 

29  o2 

20.0 

30.9 

ia.5 

20.0 
19.6 
22.6 
31.5 
21*. 0 


1 Jr 

1.175 
1.130 

1.175 

1.51 

1.585 

2.135 

2.55 

2.60 

2  o63 

2.1*6 

2.695 

1.99 
.980 
.655 
.552 
.390 
0I4O8 


T 
^ 


*S, 


60  1 

l*.l* 
60  1 

3U8o9 
3li7.9 

331.1* 
308ol 

313.9 
295.5 
308o6 
28loi* 
21*6.2 
221*. 9 
210*3 
217.6 
226o5 
215.7 


.27a 

215.2 

o2l*l 

215.8 

.213 

208o  7 

.151*5 

197.5 

TABLE  IV 

INITIAL  THEORETICAL  AND  EXPERIMENTALLY  DERIVED 
LATHI AL  STABILITY  DERIVATIVES  FCR  THE  NAVION 


T.A.S.  -     129»6  m.p.h. 

Density  Alt.    -     5000     ft. 
Gross  Wt.         «     2870     lbs* 


1 
Derivative 

Initial 
Theoretical 

Results 
of  Ref 0  5 

Final 

Experimental 
Results 

°7, 

-0.592 

-0.592 

-0.591 

% 

♦0.102 

+0.131 

+0.119 

°10 

-0.0811* 

-0.0595 

-Oo05l3 

% 

-0.1*60 

-o.ioo 

-OoU3l4 

\ 

-Oo0937 

-0.0951 

-Oo0908 

°»0 

+0,085 

+0.101 

+0.0895 

s 

-0.051 

-0.0705 

-0.106 

V 

+0.123 

+0.1U* 

+0o206 

^ 

-0.101 

-0.137 

-0,108 

C^ 

+Oo0152 

+0.0069 

+0.0128 

°n  U 

+O0OI38 

+0o0031 

+Oo00376 

Cnir 

-Oo062$ 

-0.077U 

-0.0673 

TABI£  V 

TRANSFER  FUNCTIONS  FOR  THE  NAVION  DETERMINED 
FROM  EXPERIMENTAL  FREQUENCY  RESPONSE  DATA 


T.A.S.  -     129.6  ra.poh. 

Density  Alt.   -     5000     ft. 
Gross  Wto         -     2870     lbs. 


/  x       £            -2^  s^~    20,o^!|5       -    ^.  155" 
(1)      if    = 


5*  -I-      <=\.00Z.  53  -+-    I  3.  S-4,8  S&  -4-    57.  ?&4  5  -  /•  4>3 


<5>                    -  23.470  Sz    -     /7.  z_>_TS     ~    /34>.5?? 
/p^         — ■—     ~     — 


(3) 


U»  .   24^S3  +•    M.^S^^  +  a^StS  -    22.fH 


^a  S   (S-4  +    9  _><D2  53    +-/3.5  4>?_x+57  fr<M  5    -  (.til) 


8  .  M&$3     +-    /I.9HS*    +   4/.  ?9  IS     -/.^O 

(U)        -7-    -=. 


— »— "™- 


(5)        _t>      _  3.793$*   -  n,41a"  5  3.Z.370 


<*r  3*   +    <?,60  Z  S  3    +    I3.5tf  S1  +  57.f-fV5-/.^3l 


,_.%         V>  -43  7U  S  3   -    3  3.  76  9   5  *    4    id?.  79  V  S    -  S.3~9p 

Co)         -7      ~ 
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FIGURE  7  -  MECHANISM  FOR  THE  SINUSOIDAL  OSCILLATION  OF  AILERONS  OR  RUDDER 
WITH  PROVISIONS  FOR  THE  LOCKING  OF  THE  UNDRIVEN  CONTROL 
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Figure    1 1  -  Samples   of   Raw  Data 
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Figure  12-  Samples    of  Filtered  Data 
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APPENDIX  A 

Development  of  Vector  Equations  of  Lateral  Motion 
for  Determination  of  Stability  Derivatives 


The  three  linear  equations  of  Lateral  Motion  may  be  written  in 
the  following  form: 


Side  force:  (1)       (j)  +    Cy*  ) (3     -     Cu     <p      + 


n 


Roll:  (2)     Ct*  8    +■  (D1^    CA>)  P-(j*  P*+  C*D)>=    a/rirt  0£k 
Yaw:   (3)    -C^*  (S    -  (£»    +   Cn*]>)  <t>  +  (])i+    C»*I))   V  C^Vr+  ^/j* 


The  coefficients  of  these  equations  are  defined  as  follows: 


'Y0 


70 


{    2.r  i 


cL*    -  cLfJ-r)=    9/v 


<V"  CM 


T>  P^    JLJkT  J 


CL,*-    Clr(^) 


Cn^      ■    Cn 


•" 


p  v  ^  r* 


Cnp*  "  Cnp 

On,*    -    <*, 


-2J2  r 


ZJ2T-) 


CYJ*a    CYSr    (zt) 


°\L 


Cl^r    (J5T*  ) 


cnir*-cn^  ^£L^ 


i 


X 


Clr    *  -  Ci,       /    -^-      \ 


Cn,    *-  cn  ,      /   ^- 
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The  vector  technique  is  applied  to  the  equations  in  the  following 
manner.     If  an  input      dr     to  a  linear  system  is  considered  to  have   a 
sinusoidal  variation  of  frequency   00     ,   it  can  be  considered  to  be  a 
vector  of  unit  magnitude  lying  along  the  real  axis  of  the  complex  plane 

as  shown  in  the  sketch  below. 

ImciGmary     axis 


4> 


L 


(lea  1    c  x  i  5 


>^i)24> 


60=    1.2.     radians^ 


sec, 


The  first  derivative  of  the  sinusoidal  input  is  obtained  by  multi- 
plying amplitude  by  LO  and  rotating  the  vector  through  90°.  The 
steady  state  response  of  the  system  may  be  considered  as  a  vector  of 
amplitude  ratio  R  and  having  a  direction  displaced  from  the  real  axis 
by  a  phase  angle,  (£>  .  Thus,  in  the  complex  plane,  the  vector  for 
sideslip  may  represented  as 
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In  complex  notation  the   derivative  of  either   the  control  input   or   the 
response  for  the  steady  state  condition  may  be  determined  by  multiply- 
ing by    X  u)    and  the  integral  by  dividing  by     ^,u)    « 


The  following  are  defined: 

Ay  -     \±\Co<    <|^r 
A  $  *     I   f  (Cos    $d:ir 


or 


B*  -  1  £|Sin  d'e>^r 


B  ^  -    \  J  P^S, 


B^=    )    rrl" 


$4. 


In  this  notation  with  the  input  vector  having  a  magnitude  of  unity  the 
responses  and  their  derivatives  may  be  defined  as  follows: 

(4)  $  =     A#  +  X  B* 

(5)  I) ,6   a    -oo  B/3  +  ^  to  Ar3 

(6)  DV3  -   »u>Mg--ctt)lBa 

(7)  0    -    Acpt  x  B$ 

(8)       D0     --    "(A)   B0r    ^u;A^ 

(9)  D2(D  =•  -o;2A<J)-^  u>1  b$ 

(10)      -V-     =      A  ^  +  -^  £>^i> 

(li)  D^   =   -aob-^  +U.uj  A^ 

(12)      j) z-^  =  -  uj 2-  A >l  -  >t-  <^ z  B  ^ 
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Substitution  of  equations   (4)   to   (12)   inequations    (1),    (2),   and   (3)   and 
separation  into   real  and  imaginary  parts   yields  the   following   six 
simplified  equations: 
Sideforce:    (13)    Cy6    Bfi  -   Cu    "$>$    -      -  tu(/\£  +  A^J 

(14)     Cyfl    Aa  -   C*  A«    =     Cy^r    +-  u;  (Bg  +-£>•*) 
Roll:  (15)     Q/    Bfl  +   C**«A<t>    +  ?p  u^fhj,    -   Cl?  w  Af,  =    o^fi^ 

(16)    Q*    Aj-Clju^   +    %*  u,*A>    +  Cl?(uif-^*swlA<> 
Yaw:  (17)  -CA*    B*    -  C  *  w  A<*  +  Tf*  "ZB*     *  Cfc*  ^  A^.  =    u^  B^ 

(18)  -C^    4*    4  Cn*  u:  B$  +  r^  ^A(f    -Cn*  kJ  B-*-C«*r«  WZA^ 

In  these  equations  the   A?  s  and  B   s  will  be  available  from  the 
frequency  response  curves  for  any  chosen  value  of   to    .     In  theory  it 
is  seen  that  equation  (13)   may  be  solved  for    Cyg     from  the  response 
at  a  single  frequency  since    Cj^      is  known.     Equation  (14)   may  then 
be  solved  for    Cy  ,        using  the  computed  value    Cy*    .     The  complete 
solution  of  the  yaw  and  roll  equations  may  be  accomplished  if  the 
responses  are  known  at  three  different  frequencies  since  the  maximum 
number  of  unknown  coefficients  in  any  one  equation  never  exceeds 
three.     The  method  of  averages  attempts  to  eliminate  random  error  in 
single  points  by  writing  each  equation  for  a  number  of  frequencies  and 
grouping  the  equations  thus  determined  in  as  many  groups  as  there  are 
unknowns.     The  coefficients  for  each  group  are  summed  and  the  result- 
ing simultaneous  equations  are  solved  for  the  unknowns. 

In  the  roll  and  yaw  equations  it  is  found  that  better  solutions   are 
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obtained  if  the  grouping  of  points  is  performed  in  such  a  manner   that 
all  the   points  in  a  given  group  lie  in  the  same   frequency  range   and 
further  that  the  three  groups  chosen  include  data  below,   at,   and  above 
the  dutch  roll  frequency. 

The  final  computational  equations  then  become: 


"M  Uj* 


Side  force:  (19)     Cy*     £    B^    -    C*   t  B*     =     -£    ^(Ad  +  Av) 


i^>\ 


(20)      Cv*   £~A*    -    Cu    £  A$   =   n  CY^  +     £  0,(^*6*) 


Roll:    (21a)   Q*    £  36  +    %   £«"A<|>  4   "J    £w'B^  -  C4f  £  uM^.   =    £  u>*B<t> 


(21b)    C£/6 

tC<  ^y  toe 


pie)  Q/?e*  -  e^^^^H-0^"^  ^lB* 

«^  *^>~  "^  "^  **" 

(22)     C  *T  5  /|6  -  C£  *  o>  B$  4  ^  £*M-4  +  C^f  *^i-K^  *    £^4, 
Yaw:   (23a)    -Cj"  ^B*  "  C  £"^<t  +    f*  ^  B*  +  °^   £"^Y   *    ^  3* 


Uf^r*. 


i*"*~  uj^y.  -r  tCv"~  „-■*.<:         At-      -CtolR,/ 

(23b)    -CjlU-cJf****    |#£^B#+<Ur£«At-     £    H 

(230  -  cx;  R  -  c^  £:*«+  **  P**  ♦  *?  s^*  -  i*^ 


toy*,  U>*, 


/ 
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The  development  of  the  computational  equations  was  performed 
for  the  rudder  case  with  the  coefficients    Ci  (         and    C^  ,      equal  to 
zero.     It  is  obvious  that  for  the  aileron  case  the  equations  are  identical 
to  those  developed  above,   with  the  substitution  of  Ci  <         and    C_  ( 

for    Ci  {         and    Cn  ,  and  Cv         z  0    . 

**r  ir  i  %f 

If  difficulty  in  the  solution  of  the  equations  arises  because  of  ill- 
conditioning,   it  may  be  possible  to  resolve  in  some  direction  other  than 
that  of  the  control  deflection  vector .     It  may  also  be  possible  to  elimin- 
ate some  troublesome  coefficient  between  equations  and  thus   obtain 
solutions.     No  attempt  was  made  here  to  consider  all  the  possibilities 
of  modification  of  the  equations . 

If  the  derivatives  are  known,  the  transfer  function  coefficients 
may  be  easily  developed  as  follows: 
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The  coefficients  of  the  transfer  functions  are  determined  as 
follows,    where       Cn  —       — 22:     . 


'  Ix  I 


k 


c;  =  cy*  +  a*  +  Cn*  -  a*c4  -  or*  &  -  c/e   &*> 


r«  r.i 


c»  =  c„*  t  cy;  G;  +  Cy;  c*.  +  c»*  C  -  o;  r^  -  c/r  cM; 


* 


cs'  =  c*c^  +  cx*,c4  -  o'c,*  +  cv/a.pc*-  Cvjc^c5 

-  C*   &J  C*6 

C^   -    Cv^CiJ  +  cv?pc£  -  iff  a,*  -  cr,*r    -  cWr%  ch* 


*     -l     rv*      1) 


C,     -    0,r  +    ch'r  ±* 


2 


-i.  X 


c(0'  =   c4rc/,+  cv<*  ig  cn*e  -  Cvi'c^  +  ^c* 


If 
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/  4  *  -4  4  +  -4  t  •*  *  * 


x  & 


<J  r 


^   ^  /.  .     1*  *    _±-   r*. ,      r 


Ci5     -     Cn^C**    ^  C^  Cyfl   -  CY/r0^    i£   +  CY^C 


Ti  ^Jr^** 


Ci4    ~     Cn^r  Cv^    C£p,   -   Cy^  C^  CHp  +     Cy^C^  Cv 

+"   C^h  CY6    ^>a.p 
Lis    -     C«jrCL  cM     +     Cn^  ll  O 

Once  the  transfer  function  coefficients  are  known  the  frequency 
response  may  be  determined  by  the  substitution  of     Lu)    for    s,   giving 
the  following: 
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